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We calculate the energy per particle of symmetric nuclear matter and pure neutron matter using 
the microscopic many-body Brueckner-Hartree-Fock (BHF) approach and employing the Argonne 
V18 (AV18) nucleon-nucleon (NN) potential supplemented with two different three-nucleon force 
models recently constructed to reproduce the binding energy of ®H, ®He and '^He nuclei as well as 
the neutron-deuteron doublet scattering length. We find that noire of these new three-nucleon force 
models is able to reproduce simultaneously the empirical saturation point of symmetric nuclear 
matter and the properties of three- and four-nucleon systems. 
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I. INTRODUCTION 

The important role played by the three-nucleon forces (TNFs) has been widely pointed out both in finite nuclei and 
nuclear matter calculations (see e.g. [IHl^ and references therein quoted). First indications for the inclusion of a TNF 
in the nuclear Hamiltonian arose from the discrepancy between the results of the binding energy using different 
nucleon-nucleon (NN) potentials and its experimental value. For example using high precision NN potentials, able to 
fit NN scattering data up to an energy of 350 MeV with a per datum close to 1, the ^H, ^He and ^He binding 
energies are under-predicted by about 1 and 4 MeV in the case of the three- or four-nucleon systems respectively [2l[ . 
A commonly accepted solution to this problem has been the introduction of TNF that could bridge the gap between 
the calculated binding energy based on two-body interactions and the experimental binding energies. The 

origin of such a TNF lies in the fact that nucleons are treated as point like particles disregarding therefore for their 
internal quark structure. The TNF emerges as a residual tidal force. 

In nuclear matter calculations based on non-relativistic many-body approaches similar problems arise. In such 
calculations, when only a two-body NN potential is used, symmetric nuclear matter (SMN) results over-bound and its 
empirical saturation point po = 0.16 fm“^, E/A\pg = —16 MeV cannot be reproduced. As in the case of few-nucleon 
systems, also for the nuclear matter case TNFs are considered as the missing physical effect of the whole picture. In 
addition, TNFs are likely crucial in the case of dense ^-stable nuclear matter to obtain a stiff equation of state (FOS) 
[23 - [^ compatible with the measured masses, M = 1.97 ± 0.04 Mq and M = 2.01 ± 0.04 Mq [2^ of the neutron 
stars in PSR J1614-2230 and PSR J0348-I-0432 respectively. 

In relativistic microscopic approaches, such as the Dirac-Brueckner-Hartree-Fock one, the importance of three- 
nucleon interaction is diminished [^ . In this approach TNFs are partially included by means of nucleon-antinucleon 
virtual excitations in the scalar cr-meson exchange process due to the dressed Dirac spinor in the nuclear medium. 

Although systems as nuclear matter and finite nuclei deal with the same interactions, the numerical calculations 
of the properties of these systems are commonly performed using different approaches and numerical techniques. 
Moreover, different parametrizations of the same TNF are most of the time present in literature according to the 
system treated: finite nuclei or nuclear matter. In fact, the TNF parameters can be fixed to reproduce the properties 
of few-nucleon [A = 3,4) systems, or the empirical saturation point of nuclear matter. In this paper we analyze the 
differences existing in the sector of three-nucleon interaction between these two areas of nuclear physics. We would 
like to see if the need to use different TNF parametrization is the consequence of a restricted search in the relative 
strength of some TNF terms or if it is a more fundamental problem. For example in Ref. it has been shown that 
in order to simultaneously describe the ^H, ^He and ^He binding energies and the neutron-deuteron (n-d) doublet 
scattering length it is necessary to modify some of the strengths present in the TNF. Along this line we investigate 
the possibility to find a paramatrization of the TNF suitable both for finite nuclei and many-body calculations. 

The paper is organized as follows: in the second section we present the three-nucleon force models used in this 
work; in the third section we briefly review the many-body Brueckner-Hartree-Fock (BHF) approach and we discuss 
how to include a TNF in this formalism; finally, the fourth section is devoted to show the results of our calculations 
and to outline the main conclusions of this study. 
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FIG. 1: Diagrams representing the contribution of terms a (diagram 1), h,d (diagram 2), D (diagram 3) and E (diagram 4) to 
the generic three-nucleon force W (1, 2, 3) of Eq. Q 


II. THREE NUCLEON FORCES 

The TNFs that we considered in this work are the new Tucson-Melbourne potential [3l| (hereafter TM’) and the 
three-nucleon potential based on chiral perturbation theory calculated at next-to-next-to-leading order in its local 
form [ 3 ^ (hereafter N2LOL). The TM’ potential is a revisited version of the older Tucson-Melbourne potential 
readjusted in order to satisfy the chiral symmetry. The final operatorial structure coincides with the one obtained in 
the old Brasilian three-nucleon model [Sg. These potentials, in conjunction with the Argonne V18 NN potential [3^ . 
have been recently used by the Pisa group [s^ to find a new parametrization able to reproduce simultaneously the 
binding energies of the ^He and ^He nuclei and the neutron-deuton scattering length ^Und- The TM’ and the N2LOL 
potentials can be written in the following way [s^: 


W = J2w{1,2,3) , (1) 

eye 

where IT(1, 2,3) is a generic term that can be put in the following form: 

Wil, 2,3) = aW'a(l, 2,3) + bWb{l, 2,3) + dWd{l, 2,3) + cdWd{1, 2,3) + ceWe{1, 2,3) . (2) 

In Eq. ([2]) each term corresponds to one of the different mechanism shown in Fig. [T] and has a different operatorial 
structure. The numerical values of the constants a, 6, d, ce and ce appearing in front of each term of Eq. m are 
given in Tab. U for each model. The first three terms arise from the exchange of two pions between the three nucleons 
(diagrams 1 and 2). The term a comes from ttN S-wave scattering (diagram 1) whereas the terms b and d, which 
are the most important, come from ttN P-wave scattering (diagram 2). The specific form of these three terms in 
configuration space is the following: 

IFa(l, 2, 3) = -Wo(ti • T2)(cri • r3i)((T2 • r23)y{r3i)y{r23) 

IFb(l, 2, 3) = Wo(ti • T2)[(o-i • (T2)y{r3i)y{r23) 

+ (o-i • r3i)((T2 • r23)(r3i • r23)t{r3i)t{r23) 

+ (o'! • r3i)((T2 • r3i)t{r3i)y{r23) 

+ (o-i • r23)((T2 ■ r23)y(r3i)t(r23)] (3) 

VFd(l, 2, 3) = Wo{t 3 ■ Ti X T 2 )[(cr 3 ■ 0-2 X (Ti)y(r3i)y(r23) 

+ (cTi • r3i)((T2 • T-23)(cr3 • r3i X r23)t{r3i)t{r23) 

+ (cTi • r3i)((T2 ■ r3i X (T3)t{r3i)y{r23) 

+ (<7-2 • T'23)((T3 • r23 X (Ti )y (r3i )t (r23 )] , 

Terms a, b and d are present in both the TM’ and N2L0L models. The functions y{r) and t{r) are defined below. 
Wq denotes the overall strength of these three terms and is defined in a different way in the two models. In the TM’ 
case Wq reads: 


Wo = 


Stt ruN 


( 4 ) 
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a 

b 

d 

CD 

Ce 

TM’ 

—0.87/m.n- 

-2.58/m3 

-Q.lbi/ml 

0 

0 

N2LOL 

cim^ 

C 3 I 2 

ca!'^ 

1 

-0.029 


TABLE I: Numerical values of the constants a, b, d, cd and ce appearing in front of each term of Eg. (P|l for each model. 
Constants ci = —0.00081 MeV“^, cs = —0.0032 Mev“^ and C 4 = —0.0054 MeV“^ are taken from Ref. 



b [m~^] 

d [m~^] 

Ce 

A [m,rl 

TMl’ 

-8.256 

-4.690 

1.0 

4.0 

TM2’ 

-3.870 

-3.375 

1.6 

4.8 

TM3’ 

-2.064 

-2.279 

2.0 

5.6 


TABLE II: Three different parametrizations of the TM’ three-body force. The value a = —0.87 ^ has been kept fix in all 

the three cases. See text and Ref. for details. 


where g = 14.06, and ttItt and mjv are the pion and nucleon masses, respectively. In the N2LOL model Wq is given 
by: 


W'o 




4 

2 2 


(5) 


with Ftt = 92.4 MeV and gA = 1.29. 

Term D is present only in the N2LOL model, and it provides the contribution of a two-nucleon contact term with 
the emission and absorption of a pion (diagram 3). Its local form in configuration space, derived in Ref. reads: 


ITrr(l,2,3) 


■ T 2 ) ( 0-1 • cr2)[?/(r3i)Zo(r23) 


Zo{r3i)y{r23)] 


(cTi ■ r3i)(cr2 • r3i)t(r3i)Zo(r23) 


+ (cTi • r23)(cr2 • r23)^o(?'3l)i(?'23) , 

where the constant Wq is defined as 


/ _J_\^ f rrhA.'^ f gAm^ 
{uttJ [f^J [aJ 8 


( 6 ) 

(7) 


being A^, the chiral symmetry breaking scale with a value of 700 MeV. 

Term E is also present only in the N2LOL model but not in the TM’ one. For the N2LOL model, term E gives the 
contribution of a three-nucleon contact term (diagram 4). It reads: 


where Wq is 


= W^{t3 • r2)Zo(r3i)Zo(r23) , 



( 8 ) 

(9) 



C 3 [MeV"^] 

C4 [MeV^] 

CD 

Ce 

N2LOL1 

-0.00448 

-0.001963 

-0.5 

0.100 

N2LOL2 

-0.00448 

-0.002044 

-1.0 

0.000 

N2LOL3 

-0.00480 

-0.002017 

-1.0 

-0.030 

N2LOL4 

-0.00544 

-0.004860 

-2.0 

-0.500 


TABLE III: Four different parametrizations of the N2LOL three-body force. See text and Ref. for details. The value 
Cl = —0.00081 MeV”^ has been kept fix in all the four cases. 

























4 


Although, as said before, the term E is not present in the original TM’ model, recently in Ref. this model has 
been extended by introducing a term similar to that of Eq. 

fE™'(l,2,3) = W'o^Zo(r3i)Zo(r23) , (10) 


where for simplicity the isospin dependence has been omitted. We will refer to this modification of the original TM’ 
potential also as TM’ and it is the one that we will use in all the calculations presented in this work. Note that for 
this extended TM’ force the constant ce is different from zero (see Tab. HIl) and not equal to zero (see Tab.|T|. 

The radial dependence of the five terms is encoded in the functions y(r), t(r) and For the TM’ and N2LOL 

models, the functions yir) and tir) are: 


yir) = Vo 


m 


t{r) = 


2/'(0 


( 11 ) 


where the prime symbol in /o(r) and y'{r) denotes the derivative with respect to r, the factor vo is equal to 1/3 for 
the TM’ model and to 1 for the N2LOL model. The function /o(r) is given by: 


/o(0 


127r 1 
2'k'^ 



dqq^ 


jojqr) 
q^ + ml 


FA{q) , 


( 12 ) 


with jo{qr) = sm{qr)/(qr). The cutoff function FA{q) in the TM’ model is taken as: 


while in the N2LOL model is given by: 


FA{q) 


■ 2 
- A 2 + (72 


F^iq) = exp[-q^/K‘^) 


(13) 


(14) 


A is a momentum cutoff parameter that fixes the scale of the system in momentum space. In the N2LOL, it has been 
set to A = 500 MeV, whereas in the TM’ model the ratio A/m.^. has been varied in order to describe the and ^He 
binding energies at fixed values of the constants a, b and d. In literature the TM’ potential has been used in several 
works (see e.g.. Ref. [S^) with typical values around A = 5 ttItt. 

The function Zo{r) appearing in Eqs. ®, m and OT is defined as: 


Zo{r) 


Utt 1 
m% 27r2 

TT 


dqq^jo{qr)FA{q) 


with FA{q) defined in Eq. (TT^ for the TM’ model and in Eq. (ICTl for the N2LOL one. 


(15) 


III. THE BRUECKNER-HARTREE-FOCK APPROACH 

The basic ingredient of the BHF approach in nuclear matter [11,113 is the Brueckner reaction matrix G describing 
the effective interaction between two nucleons in the presence of a surrounding medium. In the case of asymmetric 
nuclear matter [tI] with neutron density p„, proton density Pp, total nucleon density p = Pn + Pp and isospin 
asymmetry /3 = {pn — Pp)/p (asymmetry parameter), one has different G-matrices describing the nn, pp and np 
in medium effective interactions. They are obtained by solving the well known Bethe-Goldstone equation, written 
schematically as 


Gr 


1 (w) = Vr, 






-G 


Cd — €t, — 


le 


Tirj;r3T4 


(w) 


(16) 


where Tq {q = 1,2,*,/, 3,4) indicates the isospin projection of the two nucleons in the initial, intermediate and final 
states, V denotes the bare NN interaction, Qnrj is the Pauli operator that prevents the intermediate state nucleons 
(*,j) from being scattered to states below their respective Fermi momenta and w, the so-called starting energy, 
corresponds to the sum of non-relativistic energies of the interacting nucleons. The single-particle energy e,- of a 
nucleon with momentum k and mass m^. is given by 

eAk) = ^+Re[UAk)], (17) 

AtTIj- 
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where the single-particle potential Urik) represents the mean field felt by a nucleon due to its interaction with the 
other nucleons of the medium. In the BHF approximation, Ur{k) is calculated through the so-called on-energy-shell 
G-matrix, and is given by 


UAk) = Y, Y. {kkAGrr'-rr'{uj = eAk)+er'ik'))\kk')A, 
t' k'<kp ! 


(18) 



Once a self-consistent solution of Eqs. ra-dii is achieved, the energy per particle can be calculated as 

fh^k^ 


T k<kF 


A 


V 2mT 


-Re[[/,(fc)] 


(19) 


A. Inclusion of three-nucleon forces in the BHF approach 

In the microscopic BHF approach the TNFs discussed in the previous section cannot be used directly in their 
original form. This is because it would require the solution of a three-body Bethe-Faddeev equation in the nuclear 
medium and currently this is a task still far to be achieved. To avoid this problem an effective density dependent 
two-body force is built starting from the original three-body one by averaging over the coordinates (spatial, spin and 
isospin) of one of the three nucleons. The effective NN force due to the NNN one is thus 

W(l, 2) = i Tr(.^,^ 3 ) J dv 3 Y ^(1’ 2, 3) n(l, 2, 3) (20) 

eye 

In the previous expression n(l, 2, 3) is the density distribution of the nucleon 3 in relation to the nucleon 1 at ri and 
nucleon 2 at r 2 . The function n(l, 2, 3) represents the effect of the NN correlations and will suppress the contributions 
from the short-range part of IF(1, 2,3). In the following we adopt an usual choice used in literature [i^, 

n(l,2,3) = p/(l,3)5'(2,3), (21) 

where 5(1,3) and 5(2,3) are the correlation functions between the nucleons (1,3) and (2,3) respectively. The latter 
quantities can be written as 5(1,3) = 1 — 5(1,3), where 5(1,3) is the so-called defect function (and similarly for 
5 ( 2 , 3 ) ). Within the BHF approach the defect function should be calculated self-consistently with the G-matrices 
(US and the single particle potentials (fTS|) . Thus the average effective two-body force (1^ should be calculated 
self-consistently and added to the bare NN force at each iterative step of the calculations. 

To simplify the numerical calculations and following Is^l , in the present work we use central correlation functions 
g{i,j) independent on spin and isospin. Moreover, it has been shown |5ll.l52j| that this central correlation functions, in 
which are included the main contributions of the ^S'o and ^Si channels, are weakly dependent on the density, and can 
be approximated [H, by a Heaviside step function 9{rij — Vc), with rc = 0.6 fm in all the considered density 

range. Note that the average procedure has to be performed for each term involved in the cyclic permutation in Eq. 

(HD- 

In the following we report the expressions we used to perform the reduction of the original TNF to the effective 
density dependent two-body one ( lii,ll^). 

For a generic function F'(r 3 i, 7 - 23 , ■ri 2 ) where r 3 i, r' 23 , ri 2 are the lengths of the three sides of the triangle as shown 
in Fig. [2] we have: 


J *3 (Ti • f3i(T2 • fai E(r3i,r23,r-i2) cosO = 


27T 

3r 


p + 00 ^1^31—r23| 

drsi / dr23 r^i r23 E(r-3i,r23,ri2) cos9 (cri ■ CT 2 cos9 + 831 ^ 31 ) Q) 

12 Jo J 


( 22 ) 
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FIG. 2: Coordinates of three particle system. 


J dr3 a-1 • ■r3i(T2 • rsi F{r 3 i,r 23 ,ri 2 ) = 

2^ r+oo |•\r3l-r23\ 

■5- / drsi / dr 23 ^31 ^23 ^"(^31,r23,ri2) (cri 0-2 cos 9 + S 3 i{f31) Q) 

Jo 7 ^ 31 + 7-23 


(23) 


J dr3 S3i{f3i) F{r3i,r23,ri2) = 

2 ^ r+°° r|i'3i-r23| (24) 

- 512 ( 1 ^ 12 ) - / dr3i / (ir23 r3i r23 ^’(^’ 31 , ?'23, ?'12) P2{cos9y) 

Jo '^'r31+T"23 

Where cos0 = 7^31 ■ f 23 , P 2 {cos 9 ) = ^(3 cos^9 — 1), Q = cos9 — ^sin9xsin9y (see Fig. [2] for the definition of angles 
9, 9x and 9y) and S' 3 i(r 3 i) = 3 0-3 • f 3 i cri • f 3 i - 0-3 • cri. 

When we consider the term fF(l, 2, 3), the trace operator acting over 0-3 and T3 produces a factor 4 on the terms 
Wa, Wb, Wd, while makes Wd vanishing due to the traceless property of the cr matrices. On the other hand in the 
other two cyclic permutation FF(2,3,1) and VF(3,1,2) all the previous terms make zero because there is always an 
explicit linear dependence on 0-3 and T3. For the last term We we have to discuss separately the TM’ and the N2LOL 
models. The N2LOL model has a dependence on Ti ■ Tj in the term We so also in this case only the permutation 
FF(1,2,3) survives. For the TM’ model we have instead no isospin dependence in We so all the three permutations 
of (1, 2, 3) give contribution to the effective two-body force. 

Using the above formulas we can perform the two-body reduction of the original three-body force W to the effective 
two-body one. The corresponding expressions for the TM’ and the N2LOL models can be found in the appendix IVII 
The final effective two-body force is finally added to the bare NN interaction and the energy per particle is obtained 
in BHF approximation as discussed before. 

We want to stress that our average do not take into account some exchange contributions coming from closing 
a nucleonic fermion line over two different nucleons. These contributions are better evaluated starting from the 
momentum space form of the three-body potential Q. In addition another possible improvement to the average, is 
to close the fermionic line considering an interacting propagator Q. These tasks are beyond the scope of this work 
and will be considered in the future. 


IV. RESULTS AND DISCUSSIONS 

We now present the results of our calculations of the energy per particle of symmetric nuclear matter (SNM) and 
pure neutron matter (PNM) using the AV18 NN potential supplemented with the TM’ or N2LOL three-nucleon force. 
Making the usual angular average of the Pauli operator and of the energy denominator [i^ . , the Bethe-Goldstone 

equation m can be expanded in partial waves. In all the calculations performed in this work, we have considered 
partial wave contributions up to a total two-body angular momentum Jmax = 9. 
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FIG. 3: (Color on line) Energy per particle EjA of symmetric nuclear matter as a function of the nucleonic number density 
p for the three parametrizations of the TM’ model (left panel) and the four parametrizations of the N2LOL one (right panel) 
considered in the present work. The green double-dash-dotted line, in both panels, represents the energy per particle with 
no three-body force contribution and using the AV18 NN potential. The empirical saturation point of nuclear matter po = 
0.16 ± 0.01 fm“®, E/A\pg = —16.0 ± 1.0 MeV is denoted by the green box in both panels. 


Following Ref. [s^, we consider the three parametrizations for the TM’ model reported in Tab. HIl (hereafter called 
TMl’, TM2’ and TM3’), and the fonr parametrizations of the N2LOL model reported in Tab. IIIII (hereafter called 
N2LOL1, N2LOL2, N2LOL3 and N2LOL4). 

In Fig. [3] we show the energy per particle E/A of symmetric nuclear matter. The green double-dash-dotted line, 
in both panels, represents the energy per particle with no three-body force contribution. The resulting saturation 
point is Po = 0.23 fm“^, E/A\q = —16.43 MeV, to be compared with the empirical saturation point of nuclear matter 
Po = 0.16 ± 0.01 fm“^, ElA\„r. = —16.0 ± 1.0 MeV (green box in both panels of Fig. [S]). We next introduce the 
three-body forces of Ref. using the average procedure described in the previous section. In the case of the TM’ 
model (left panel) the three-body force produces a sizeable repulsive effect (be. E/A increases with respect to the case 
with no TNF) in all the considered density range, and shifts the calculated saturation point (see Tab. IIVI) to a density 
lower than the empirical one. At the empirical saturation density E/A increases by /S.E = 4.9 MeV (6.8 MeV) in the 
case of the TMl’ (TM3’) interaction. At twice the empirical saturation density, be. p = 0.32 fm~^, /S.E = 20.9 MeV 
(28.1 MeV) in the case of the TMl’ (TM3’) interaction. 

The outcome is notably different in the case of the N2LOL (right panel) three-body forces. In this case the TNF 
produces a decrease of E/A in all the considered density range. At the empirical saturation density (at twice the 
empirical saturation density) E/A decreases by AA = —1.4 MeV (—2.9 MeV) in the case of the N2LOL1 interaction. 
The contrasting effect on the energy per particle of SMN of the two TNF models, illustrated in the two panels of Fig. 
|3l is mainly due to the different action of the repulsive component on the two-body effective force W(l,2) derived 
from the genuine TNF. As discussed before, the N2LOL model has a nontrivial isospin dependence in the repulsive 
term We (see Eq. (1551) ). If the parameter Ce is positive (negative) the final contribution is repulsive (attractive) on 
channels with isospin T = 1 but is attractive (repulsive) on channels with isospin T = 0. On the other hand in the 
TM’ model there is no isospin dependence on the repulsive part of the three-body force (see Eq. (|2^ ) so We gives in 
all channels a repulsive contribution. 

The value of the saturation density, energy per particle and symmetry energy at the saturation density are reported 
in Tabs. 113 and m for the TM’ and the N2LOL models respectively. Both models fail to reproduce the empirical 
saturation point of SNM. This is not surprising (see e.g. Hi) since in the present nuclear matter calculations we used 
TNF models whose parameters have been determined to reproduce the properties of light (A = 3, 4) nuclei and 
the neutron-deuteron doublet scattering lenght. 
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po (fm”’^) 

E/A (MeV) 

Esym (MeV) 

no TNF 

0.23 

-16.43 

35.25 

TMl’ 

0.12 

-11.11 

26.55 

TM2’ 

0.11 

-9.98 

22.49 

TM3’ 

0.11 

-9.96 

23.05 


TABLE IV: Saturation properties of symmetric nuclear matter for three different parametrizations (first column) of the TM’ 
three-body force. The entry ”no TNF” refers to a calculation without three-body force and using the AV18 NN potential. The 
other entries in the table are: the saturation density (second column), the value of energy per particle at saturation (third 
column) and the value of the symmetry energy at saturation (forth column). 



Po (fm“®) 

E/A (MeV) 

Esym (MeV) 

N2LOL1 

0.26 

-18.47 

42.30 

N2LOL2 

0.25 

-17.90 

40.02 

N2LOL3 

0.26 

-18.09 

41.06 

N2LOL4 

0.25 

-16.90 

36.25 


TABLE V: Saturation properties of symmetric nuclear matter for three different parametrizations of the N2LOL three-body 
force. 


In Fig. S] we plot our results for the energy per particle of pure neutron matter. The green double-dash-dotted 
line, in both panels, represents EjA when the TNF is not included. In the case of the TM’ model (left panel) 
the TNF produces a sizeable repulsive effect, in all the considered density range, as compared to the case with no 
TNF. For example, at the empirical saturation density (at twice the empirical saturation density) EjA increases 
by = 6.7 MeV (25.5 MeV) in the case of the TMl’ interaction. The effect of TNF on the energy per particle 
of PNM is less pronounced in the case of the N2LOL model (right panel). In the particular case of the N2LOL4 
parametrization, TNFs make PNM softer with respect to the case where TNFs are not included. Notice that in the 
case of pure neutron matter, we have just the contribution of the T = 1 isospin channel and therefore, also for the 
N2LOL model one has a pure term that provides repulsion. Nonetheless the strength of term associated to repulsion 
in the TM’ model is stronger than the corresponding one for the N2LOL and consequently a stiffer neutron matter 
equation of state is obtained. 

The nuclear symmetry energy is dehned by 


^sym {p) 


1 d^E/A 

2 5/32 



(25) 


where E/A, the energy per particle of asymmetric nuclear matter, is calculated using Eq. (1191) within the BHF 
approximation. 

It has been numerically demonstrated by the authors of Ref. jssj and afterwards confirmed by various microscopic 
calculations (sgI - I^ adopting different realistic NN interactions, that the energy per particle of asymmetric nuclear 
matter can be accurately reproduced by the following relation: 


(P) ~ A Esym{p)P 


(26) 


Thus, in good approximation, the symmetry energy can be expressed in terms of the difference of the energy per 
particle between symmetric {j3 = 0) and pure neutron matter (/3 = 1): 


i^.ym(p) = |(p,l)-|(p,0). (27) 

The symmetry energy for the TM’ and the N2LOL models is shown as function of the density on the left and right 
panels of Fig. [5] respectively. The green double-dash-dotted line, in both panels, represents the symmetry energy with 
no three-body force contribution and using the AV18 NN potential. In Tabs. HVl and IVl we have reported the values 
of the symmetry energy for the two TNF models at their respective calculated saturation points po (second column 
in Tabs. m and 0. 
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FIG. 4: Energy per particle EjA of pure neutron matter as a function of the nucleonic number density p for the three 
parametrizations of the TM’ model (left panel) and the four parametrizations of the N2LOL one (right panel) considered in 
the present work. The green double-dash-dotted line, in both panels, represents the energy per particle with no three-body 
force contribution and using the AV18 NN potential. 



Esym (MeV) 

L (MeV) 

no TNF 

28.79 

51.3 

TMl’ 

30.14 

53.6 

TM2’ 

29.30 

50.3 

TM3’ 

28.65 

48.4 

N2LOL1 

31.14 

56.6 

N2LOL2 

30.59 

53.7 

N2LOL3 

30.65 

54.1 

N2LOL4 

28.92 

44.8 


TABLE VI: Symmetry energy and slope parameter L at the empirical saturation density pnm = 0.16 fm ^ for different TNF 
models. 


To compare our results with the value of the symmetry energy extracted from various nuclear experimental data 
[HH, we report in Tab. ED Esym calculated at the empirical saturation density po = 0-16 fm ^ for the different 
TNF models considered in this work. In the same table, we also report the so called slope parameter 




dp 


PQ 


(28) 


As we can see (Tab. IVII) our calculated Esym and L lies within the ranges of values extracted from experimental data 
@|: Esymipo) = 29.0 - 32.7 MeV, and L = 40.5 - 61.9 MeV. 

The TNFs of Ref. have been recently employed for SNM and PNM calculations in using both the variational 
and the auxiliary field diffusion Monte Carlo approaches. The authors of Ref. used the Argonne V8’ NN 
potential which is a simplified version of the AV18 potential [s^ . truncated after the linear spin-orbit components 
and refitted to have the same isoscalar part of the AV18 in all the S and P waves, as well as in the ^Di wave and 
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FIG. 5: Symmetry energy as a function of the nucleonic number density p for the three parametrizations of the TM’ model 
(left panel) and the four parametrizations of the N2LOL one (right panel) considered in the present work. The green double¬ 
dash-dotted line, in both panels, represents the symmetry energy with no three-body force contribution and using the AV18 
NN potential. 


its coupling to the wave. Notice that also in no one of the TNF models was able to reproduce the correct 
saturation point of symmetric nuclear matter. A direct comparison of our results with those reported in would 
be ambiguous and inconclusive for the following reasons. First of all due to the difference in the two-body used in the 
calculations: the AV18 in the present paper and the AV8’ in [b^. Second, as discussed in (^,11^, the use of different 
many-body approaches affects considerably the results particularly in the case of SNM. 


V. CONCLUSIONS 

A new generation of TNF models has been recently proposed in Ref. [s^. These new TNFs have been used, in con- 
juction with the Argoime V18 two-nucleon interaction, and their parameters have been determined to simultaneously 
reproduce the measured binding energies of ^H, ^He and ^He nuclei as well as the measured n-d doublet scattering 
length. A correct prediction for these physical quantities can be regarded as a severe requisite for a realistic nuclear 
Hamiltonian containing two- and three-nucleon interactions. As shown in Ref. [s^, this requirement was not fulfilled 
by several of the TNF models available in the literature. 

In the present work, we have calculated the energy per nucleon of symmetric nuclear matter and pure neutron matter 
within the BHF approach and using the same nuclear Hamiltonian as the one used in Ref. {i.e. without changing 
the original values of the TNF parameters) with the purpose to test this Hamiltonian in a many-body context. We 
found that none of the TNF models given in [30l| is able to reproduce the empirical saturation point of symmetric 
nuclear matter. This outcome concords with the results obtained with the Urbana VH [b^ TNF when used in few- 
body nuclei and nuclear matter (see also [3). In particular, in the case of the AV18-I-TM’ Hamiltonian, both the 
calculated saturation density and the corresponding binding energy per nucleon {B/A = —EjA) are understimated. 
The TM’ model for the TNF produces a strong repulsive effect both in SNM and PNM, in all the considered density 
range. In the case of the AV18-I-N2LOL Hamiltonian, the TNF produces a decrease of E/A (an increase of the binding 
energy per nucleon) in all the considered density range. 

The reasons why TNFs fitting few-nucleon systems are not able to reproduce the empirical saturation point of 
SNM can be various. First of all in the BHF approach we are still not able to use the genuine three-body force 
but, due to the technical reasons explained above, we are forced to include the TNF contribution by performing the 
average procedure to generate an effective density dependent two-body force. In this way some terms get lost in 











11 


the average and therefore, in a way, we are not dealing with exactly the same three-body force employed in finite 
nuclei calculations. An interesting test could be the use of our effective two-nucleon force in finite nuclei calculations. 
So doing, we may understand at a deeper level, the differences between the exact procedure in which we employ 
the genuine three-nucleon force and our simplified approach that involve an average force. We plan to perform this 
interesting comparison in a future work. In addition, the inclusion of the exchange terms in the average or the use of a 
correlation function explicitly dependent from spin and isospin will be definitely an improvement of our calculations. 
These aspects will be considered in a future work. 


VI. APPENDIX 


The explicit expressions for the TM’ and the N2LOL models are reported in the following. For the TM’ model we 
have: 

iri2) = I ("^1 ■ "^ 2 ) [o-i • (T2 Vairi2) + Si2iri2) Vtiri2) + Wr(D2)] ■ (29) 


27r /■+“ rl’'3l-r-23l 

Va{ri2) = - / drai / dr 23 r^i r23 [ - 5 —?'3i ?'23 2/(?'3i) y{r23) cos6+ 

^12 Jq J r3i-|-r23 


b Wo r 3 i r23 y(r3i) y{r23) + rfg cos^O + ^(rgi) y{r23) + rfg t{r23)y{r3i) ) 


b W( 


(30) 


9^(r3i) g^(r23) ■ 


27r f+°° /■ 1^31-1-231 Q 

Vt{ri2) = - / drsi / dr23 rsi r23[ -^»'3i »"23 2/(?'3i) y{r23) Q+ 

‘f’12 Jo Jr3i-\-r23 ^ 

^(»' 3 i) t{r23) cos9 Q + P 2 {cosey) rl^ t(r 3 i) y{r23)+ 

P 2 {cos 0 ^) rig t{r23) y{r3i) )] 3^(?’3i) 5^(?’23) • 

2 -j^ r+°° r\r3l-r23\ 

Vr{ri2) = Ce Wo — / droi / dr23 rsi r23[ 2;o(?'3i) ^o(?'23)+ 

^12 Jo J r3i+r23 

Zoiroi) Zo{ri 2 ) +zoiri 2 ) Zoijr 23 ) ] g'^{r 31 ) g^{r 23 ) ■ 

For the N2LOL model we have: 


(31) 


(32) 


1^^2^0^(ri2) = ^ (D ■ T 2 ) [cTi • CT2 Va(ri2) + £' 12 (^ 12 ) f’t(D2) + Vr{ri2)] ■ (33) 


27r /■|'■3l-r■23l 

Va{ri 2 ) = - / dr 3 i / dr23 rsi r23 [ - 5 —?'3i ?'23 2/(?’3i) y{r23) cos9+ 

f ’12 Jo dr3i+r23 

b Wo r 3 i r23 y{r3i) y{r23) + ^-^^31 ^23 -t ^^(r|i t(r 3 i) y{r23) + r^g t(r23)y{r3i))+ 
( r|g t(r3i) 2:(r23) -Hr|g t{r23)z{r3i) ) ] g^rgi) g^(r23) ■ 


27r /■+°° r\r3i-r23\ ^ 

Vt{ri2) = - / drgi / dr23 rgi r23[ - f ’31 f ’23 y{r3i) y{r23) Q+ 

ri 2 Jo Jt3 


/o 
h VFo 


/r3i+r23 
^2 \ 2 


rii t{r3i) rig t(r23) cos9 Q + P2{cos9y) r|g t(r 3 i) y(r23)-|- 


( 35 ) 


P2{cos9^) r|g t(r23) y(r 3 i) )-h 


( ^31 ^(^ 31 ) z{r23)P2{cos9y) + r 23 t{r23)z{r3i)P2{cos9^) ) ] g (rgi) g (r23) ■ 


3 
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i 


|r-3l-'r23l 


dr23 rsi r23[ Zo{r 3 i) Zo(?'23) ] g'^{r 31 ) g'^{r 23 ) ■ (36) 
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